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Equation by Finite Differences' 

Wolfgang Wasow 

The difference between the solution of the Laplace differential equation and the Laplace 
difference equation, defined in the same rectangle and assuming the same boundary values, 
is estimated under the assumption that the boundary function possesses a bounded third 
derivative. The bound obtained is of the order of magnitude of the square of the mesh length. 



1. Introduction 

Probably the most widely used approach for the 
numerical solution of differential equations for which 
simple analytic solutions are not available consists in 
replacing the differential equation by a finite differ- 
ence equation. The difference between the exact 
solutions of these two problems is usually called— 
somewhat inappropriately — the truncation error. 

It is highly desirable, for numerical calculations, 
to have adequate estimates of this truncation error. 
There exist in the literature numerous investigations 
of this problem. We shall be primarily concerned 
with the case of the first boundary value problem for 
elliptic partial linear differential equations. Some 
references to the literature for this case are listed at 
the end of this paper. All these investigations, as 
far as they are mathematically complete, suffer from 
the defect that the estimates given depend on bounds 
for the derivatives of the unknown solution of the 
differential equation itself. Such bounds can in gen- 
eral not be found without effectively solving the dif- 
ferential equation. Thus the value of these estimates 
consists essentially in giving information about the 
order of magnitude of the truncation error in terms 
of the mesh length used, and in helping to form 
reasonable guesses as to the size of the error. 

It turns out that, with the difference expressions 
used in most numerical work, and for a square net of 
mesh length h the truncation error is of the order 
0(h^), provided the boundary values and the equa- 
tions of the boundary curve C possess bounded Jourth 
derivatives, and provided the boundary values of the 
difference equation problem are suitably adjusted 
(see, e. g. [31^). This estimate is then vahd in the 
whole domain B. 

The bound for the truncation error given by P. C. 
Rosenbloom in [6] is somewhat different. It depends 
only on the modulus of continuity in B+C of the 
unknown harmonic function, but it is a pointwise 
estimate that deteriorates indefinitely at points near 
the boundary. 

Thus there remain a number of unsolved problems 
of theoretical as well as computational interest: The 
ideal would be a pointwise estimate of the truncation 
error, valid if the boundary curve and the boundary 
values are piecewise analytic. This error estimate 
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should have the highest possible order in h, and, at 
the same time, it should not be unduly wasteful. It 
should, of course, depend only on quantities that can 
be easily computed from the data. 

From the theoretical viewpoint the order of mag- 
nitude of the error is the most interesting quantity. 
In particular, one would like to know the exact way 
in which a weakening of the smoothness required of 
the data affects the global and local order of the error. 

For computational purposes an estimate will be 
most useful if it yields small bounds for the error 
at moderate sizes of h. Thus lOh is a better bound 
in practice than 10,000^1 

The contents of this note are intended as a first 
exploratory step toward the examination of these 
questions. The error estimate for Laplace^s equa- 
tion in a rectangle given below in formulas (18) and 
(19) is of order 0{h^), with coefficients of moderate 
size, although it depends on bounds for the third 
derivatives of the boundary values only, and in spite 
of the fact that a rectangle is a domain with corners. 
It is hoped that this special result may help to 
stimidate more general investigations. 

2. Truncation Error for Dirichlet's Problem 
in a Rectangle 

Let B be the rectangle with vertices (0,0), (1,0), 
(1,6), (0,6) where 6 is a rational immber. The 
prescribed boundary function / is assumed to be 
continuous on the boundary C of B and to possess 
bounded third derivatives on each closed side of C. 

We denote the Laplacian, as usual, by Au and 
write AnV for the expression 

p[?;(x + /t,2/)+^(2^-/^,2/) + «^(a:,]/+A') 

+ v(x,y — h)-4.v(x,y)]. 

For any value of h such that 1/h and b/h are integers 
we denote by Uh(x,y) the function, defined at all 
netpoints (nh, mh) in B-\-C, where n,m are integers, 
for which 

Aj{iin=^, in B 
and 

Uh=f, on C. 

If u designates the solution of the corresponding 
Dirichlet problem for the operator A, then we are 
concerned with the truncation error u—Un. 



345 



It is desirable to deal with boundary values that 
vanish at the four vertices. To this end we intro- 
duce the harmonic polynominal Q{x,y), defined by 
h'Q{x, y)=A^o (x—1) {y—b)—AioX {y—h)+A^i,xy— 
AQb{x—l)y,^ where ^oo, Aio, etc. are the values of/ 
at the vertices. This function solves both the differ- 
ence and the differential equation. If we replace the 
boundary values / by 

we obtain therefore a new problem for which the 
truncation error is the same as for the original prob- 
lem, since the truncation error corresponding to the 
boundary problem determined by the values of 
Q{x,y) on C is zero. It is clear that/* vanishes at 
the vertices. Also the second and higher derivatives 
of/* are the same as those of/. 

The truncation error can be considered the sum of 
four terms, each corresponding to boundary values 
that are zero except on one of the four sides of the 
rectangle. Therefore, we assume temporarily that 
/* is zero except on the side y=0. There /*=/* (x) 
is a function of x alone. 

The solutions of Au=0 and AnUn^O for these 
boundary values can be written, respectively, 



u{x,y) = ^Cng(y,nT) sin nwx, 



l/h 



Un{x,y) = ^ yng(y,fin/h) sin mrx, 



where 



^, = 2jJ/*(0si 



sin mrtdt 



l/h 

Tn==2/iXl/*(^'^) sin JiTrh 

r=l 

g(y^z)=sm\i {b — y)2lsm\i hz, 
and jSri is the solution of 

smh ^=sm —^ — 



(1) 

(2) 

(3) 

(4) 
(5) 

(6) 



These well-known formulas can be readily verified 
by substitution. 

We next collect a few facts needed for the com- 
parison of u and u^. 

To find an estimate of the Fourier coefficients c^ we 
integrate the right member of formula (3) by parts 
three times, obtaining the formula 



{ni^f 



[(-l)''-V"(l)-/"(0)] 
, 2 



This implies 



|c„|<X«,-^ 



/>"») 



cos nirtdt. 



(7) 
(7) 



where 

K=-, |[/"(1)1 + I/"(0)|] + A max |/'"(x)|. (8) 

TT X 

The Fourier coefficients Cn and the ' 'interpolation 
coefficients^', jn, are related by the equation 



Jn — (^n'ir2u \(^2k/h+n — (^2k/h-n)j (^ — 1> 2, 



■,llh). 



This relation can readily be proved by expressing 
f^irh) in (4) by its Fourier series. (See also the 
references in [7].) Combining (7) with the last for- 
mula we find, iov n<l/h, 

\Cn-yn\<K^ {(2k/h+n)-' + {2klh-n)-'] 



Hence, 



\Cn-yn\<l^7Kh^, forn<l/h. 



(9) 



Next we need an estimate of the difference g(y, nw) — 
9 (y, ^n/h). To this end we express pn/2, as given 
by (6), in terms of nh7r/2 by means of Taylor's for- 
mula. A straightforward calculation shows that 



/nhirV 



R 



~2~~Y" 

where R is the value of 

-| cos <T (1-2 sin^ (t) (1 +sin2 (j)-^i^ 

o 

at a certain mean value of a in the interval 
0<o-<?^/i7r/2. 

\fin — nhTr\<{nhTryil2. 



Hence 



(10) 



In order to estimate the difference g{y,n7r) — 
g {y, finlh) with the help of the theorem of the mean 
we have to estimate the right member of the formula 

X[{2h — y)^m\iyz — y^m\i{2h — y)z]. (11) 

We observe that sinh tjt is an increasing function of 
t, for ^>0, Therefore, 



d 



'V, 



<-T- g{y, z) <| sinh-2 hz sinh {2b— y) z, 

for0>O, ?/>0. (12) 
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To this expression we shall apply the inequalities 

sinh A:^<6'^*-^)'sinh^, forO<A:<l, (13) 



and 



smht>^il-e-''^)e\ ^>fi>0. (14) 



The second of these inequalities is obvious. The 
first follows from 

Also, we shall need the inequalities 



1^ 'n.<^nlh<Trn, forl<n<l/A. 



(15) 



The second of these is a consequence of (6); the first 
can be proved as follows: 



Pn "nJiTT / . nJiT Inhw 



and 



/ . nliT Inh7r\. 



nhw 2 

T TT 



sinh^==-T^ ( sin— ^r-/ -^T- )>-7r- -=nh. 



Since sinh xjx is an increasing function of x for x = 0, 



Hence 



sinh^/f^<f/arsinh|<|. 



|>|sinh|>|nA, 



which completes the proof of (15). Applying (13), 
(14) and (15) to (12) we obtain 



d 



l<n<llh 



.0<y<b 
and therefore, using (10), 

g(y> y- g(y, ^^)| < (i-e-''^r' jj ye-^-^Vh^ 

torl<n<l/h. (16) 

We conclude these preparatory remarks with the 
inequality 



0<g(y,z)<l, ioTO<y<b, 



(17) 



which is an immediate consequence of the definition 

of g{y,^)^ 

We are now ready to estimate \u—Uh\. To this 
end we write, using (1) and (2), 



\u-u,\<R,+R2+Rs, 



1//1 



^i=S \(^n—yn\g{y,mr) 



l/h 
n=l 

^3= X) \(^n\giy,nT). 

n=l-\-l/h 

From (9) and (17) we conclude that 

Ri<1.7Kh' 
and from (7) it follows that 

In order to estimate i?2, we note that, in view of (7) 

and (9), 

\yn\<K(n-'+1.7h'). 
Using also (16) we obtain, by summation, 

R2<(l-e-''^')-'K(2.7)^yie''^'-l)-'h'. 
1 z 



Since 



yKe' 



-l)<j^ ioTy>0, 



this yields 

R2<3.24:{l-e-''^')-'Kh' 
and, hence, adding these inequalities. 

If we denote by AI2 the numerical maximum of the 
second derivative of the given boundary function at 
the vertices of the rectangle, and by M3 the numeri- 
cal maximum, anywhere on C, of the third derivative 
of the boundary function, the last inequality can be 
replaced, in view of (8), by 

!ti-'?/,i<[.297 + .676(l-6-«^/^)-i][M2 + .319M3]/il 

Finally, if the sides of the rectangle are A, B in- 
stead of 1,6 a homothetic transformation yields the 
formula 

|^-^/.|<[.297 + .676(l-6-«^/^^)-i][M2 + .319.4M3]^l 

We recall that this formula is valid only for bound- 
ary values that vanish except on the side ^=0. For 
the full truncation error three analogous expressions 
have to be added to the right member of the last 
inequality, leading — with some rounding off to 
shorten the expression — to 

\n-Un\<l^ + lA{l-e-^^'^^)-'][M2+.Z2AM,W 
+ [.6 + 1.4(1 -6-8^/3^-1] [M2+.325M3]fi,2. 

(18) 
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For the square of side A this reduces to 

\u-Un\<(4:.3M2+lAAMs)h'. (19) 

Added in proof: The function 

<(^,2/) = X) CnOiy, 13 n/h) sin nirx 

n = l 

coincides with Un (x, y) at all net points. Its use in- 
stead of Ufi {x, y) simplifies the calculations somewhat 
and leads to better estimates for the truncation error, 
namely 

\u-ut\ < [.13+. 5 (1 -e-8^/3^)-'] [M2+0.32^M3] h' 



and 



+ [.13 + .5(l-6-«^/^^)"'][M2+0.325M3]A2 
\u-ut\ <{lAM2 + 0A?>AM^)h^ 



instead of (18) and (19). 
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